Transient solar oscillations driven by primordial black holes 
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Stars are transparent to the passage of primordial black holes (PBHs) and serve as seismic de- 
tectors for such objects. The gravitational field of a PBH squeezes a star and causes it to ring 
acoustically. We calculate the seismic signature of a PBH passing through the Sun. The background 
for this signal is the observed spectrum of solar oscillations excited by supersonic turbulence. We 
predict that PBHs more massive than 10 21 g (comparable in mass to an asteroid) are detectable 
by existing solar observatories. The oscillations excited by PBHs peak at large scales and high 
frequencies, making them potentially detectable in other stars. The discovery of PBHs would have 
profound implications for cosmology and high-energy physics. 



The existence of black holes is one of the most startling 
predictions of general relativity. Astronomers have dis- 
covered two populations of black holes: stellar-mass 
(n^BH ~ 1OM ) black holes that form in the collapse 
of massive stars and supermassive black holes (10 6 M Q < 
mBH ^$ 10 10 M Q ) that reside in galactic centers. How- 
ever, general relativity allows black holes to have any 
mass. Black holes much less massive than a solar mass 
Mq could have formed from density perturbations in the 
early Universe. Such perturbations were created with 
a wide range of wavelengths and amplitudes. Galaxies 
slowly grew from perturbations whose amplitudes were 
initially very small. Density perturbations with higher 
initial amplitudes might have gravitationally collapsed 
into primordial black holes (PBHs) [1 . The mass m bh of 
such PBHs reflects the mass contained within the parti- 
cle horizon of the Universe at the time they were formed. 
PBHs as small as the Planck mass mpi = yJhc/G ~ 1CT 5 
g may have formed, but those with masses less than 
Wevap — 5 x 10 g would evaporate in less than the 
age of the Universe [2]. 

A density perturbation will collapse into a black hole if 
its self-gravity exceeds its pressure support [51 0] . When 
this pressure is supplied by radiation as in the early Uni- 
verse, PBHs of any mass are equally likely to form if there 
is a flat power spectrum of primordial density perturba- 
tions (spectral index n s ~ 1) as indicated by observations 
of the cosmic microwave background [5J- PBH produc- 
tion may be greatly enhanced at a particular mass scale 
if the pressure were suddenly reduced, such as during the 
QCD phase transition [BJ. The discovery of PBHs of a 
given mass would thus provide insight into high-energy 
physics at the temperature at which this mass was con- 
tained within the particle horizon. 

PBHs are also of great interest to cosmology. They are 
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collisionless and nonrelativistic, making them ideal dark- 
matter candidates. Observational constraints on the cos- 
mological density £Ibh of PBHs depend on their mass 
m BH- PBHs with masses slightly above the evapora- 
tion limit m eV ap emit Hawking radiation, including 7-rays 
with a spectrum peaking around 100 MeV [7]. Observa- 
tions of the extragalactic 7-ray background by the En- 
ergetic Gamma Ray Experiment Telescope (EGRET) [8 
set an upper limit flsn < 5 x 10~ 10 for rngfl = "Vvap 
[9] . Since the luminosity and temperature of a PBH scale 
as ^^Rf7 

and m BH - respectively, it is more difficult to 
observe Hawking radiation from larger PBHs. PBHs of 
mass mBH ^ 10 17 g could constitute the entirety of the 
dark matter {Qbh = &dm — 0.23) without violating 
observational constraints on Hawking radiation. 

Microlensing surveys constrain the abundance of more 
massive PBHs. If a PBH passes between an observer 
and a background star, that star will be gravitationally 
microlensed, briefly increasing its observed flux. The du- 
ration of each lensing event is proportional to the square 
root of the lens mass, implying that a survey with finite 
cadence will miss lenses that are too small. The EROS 
(Experience de Recherche d'Objets Sombres) microlens- 
ing survey sets an upper bound of 8% on the fraction of 
the Galactic halo mass in the form of PBHs with masses 
in the range 0.6 x 1O" 7 M <m BH < 15M [TO! [IT]. 

Collisions with Galactic stars might constrain PBHs 
in this permitted range 10 17 g < nigs 10 26 g |12j . 
PBHs passing through a star deposit energy by dy- 
namical friction. The fractional loss of kinetic energy 
/ - 10~ 7 (m Bff /10 26 g) is tiny for PBHs with speeds 
comparable to the velocity dispersion of Galactic dark 
matter. PBHs will therefore pass through stars without 
slowing down or accreting appreciably. A 10 20 g PBH 
was estimated to emit 1 keV x-rays with a luminosity 
of a 10 22 erg/s when passing through a main-sequence 
star like our Sun [T^] , but most of these x-rays will ther- 
malize and this luminosity is much less than the Sun's 
background x-ray luminosity L x > 10 26 erg/s [13] . 

We instead propose searching for PBHs by the distinc- 
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FIG. 1: Radial velocities V BH induced by a KT 10 M Q PBH 
passing through the Sun along the z axis. The gray sphere 
at R < 0.24_R Q is excluded from our simulations. The colors 
indicate the radial velocity of the fluid weighted by the square 
root of the density. Wave amplitudes are linearly proportional 
to the black-hole mass, so we may rescale our simulations to 
any desired value of ttibh- 



tive oscillations they excite when passing through the 
Sun. This signal was considered previously but not cal- 
culated explicitly [T2]. We simulate the generation and 
propagation of acoustic waves in the Sun by solving the 
linearized Euler equations [H] in a spherical shell with 
the moving PBH acting as the source. We solve these 
equations on a three-dimensional grid consisting of 1024 
longitudinal, 512 latitudinal, and 425 radial points ex- 
tending from t/Rq = 0.24 to 1.002 [13 [T£]. Our sim- 
ulations exclude the core (r < 0.24 R®) because of the 
coordinate singularity at r = 0. We treat the PBH as a 
ball with radius A comparable to our grid spacing. Al- 
though we cannot resolve the PBH's Schwarzschild ra- 
dius Rs — 1.5 x 10 _7 (msij/10 21 g) cm, we verify con- 
vergence by repeating simulations with different values of 
A. Wavefield velocities are extracted 200 km above the 
photosphere (r — Rq) to mimic observations. 

A snapshot of one of our simulations is shown in Fig. [T] 
We can place the PBH's orbit (shown by gray dots) in 
the x — z plane without loss of generality because the Sun 
is nearly spherically symmetric. PBH orbits in this plane 
are fully specified by two parameters: the energy E and 
angular momentum L per unit PBH mass. The PBH or- 
bit in the simulation shown in Fig.[l]is parabolic (E = 0) 
and radial (L = 0). We consider only unbound orbits 
(E > 0) as gravitationally bound PBHs are extremely 
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FIG. 2: Power spectra Pi(v) of our simulated PBH signal 
V BH (left panel) and the turbulently driven noise V" ob- 
served by HMI (right panel) as functions of frequency v and 
harmonic degree I. Each discrete ridge in both power spec- 
tra corresponds to a set of p modes with the same number 
of radial nodes. The signal peaks near I ~ 40, v ~ 5.5 mHz, 
at larger scales and higher frequencies than the noise back- 
ground, helping us to distinguish them. 

unlikely. The PBH begins at Zi ~ 3Rq with an in- 
ward velocity of v z — — \J GMq Jz~i, falls radially inwards 
through the Sun's center, and ends at Zf ~ — 10i? Q . The 
total elapsed time is 8 h, with the snapshot in Fig. [T] be- 
ing taken 7.9 h into the simulation. The radial velocity 
V BH as a function of time t and angular position (8, cf>) on 
the Sun's surface constitutes the signal for our proposed 
PBH search. 

Photospheric velocities are observed by measuring the 
Doppler shift of solar absorption lines. The Helioseis- 
mic and Magnetic Imager (HMI) [T7] onboard the Solar 
Dynamics Observatory currently performs such observa- 
tions. The Sun has a discrete spectrum of global acoustic 
oscillations known as p modes because pressure provides 
the restoring force [14,. P modes driven by near-surface 
supersonic turbulence provide the dominant contribution 
to observed photospheric velocities, and constitute the 
noise background for our PBH search. 

P modes can be clearly identified after the photo- 
spheric velocity field has been Fourier transformed into 
a function of frequency v and decomposed into spherical 
harmonics Yi m {6, 4>). We show the power spectrum 

ra— — l 

of our simulated signal V BH and the observed noise V n 
in Fig. [2] The power spectrum in the left panel is pro- 
duced from the simulation shown in Fig. [l] while the 
noise power spectrum in the right panel was prepared by 
using 8 h of publicly available HMI data [HJ. The dis- 
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FIG. 3: Orbits of the PBHs in our simulations. The green 
vertical line along the z-axis shows the two radial orbits (RFF 
and RGC). The other two orbits are parabolic (E = 0). The 
inner black curve shows the orbit which spends the maximum 
amount of time within the Sun (MT), while the outer tan 
curve shows the orbit which skims the surface (SS). 



Crete ridges seen in both panels correspond to p modes 
with the same number of radial nodes. The power spec- 
trum of oscillations excited by the PBH peaks at higher 
frequencies and lower harmonic degree I than the noise 
background driven by the Sun's natural turbulence. 

This dissimilarity between the spatiotemporal depen- 
dence of the signal and noise helps us detect small sig- 
nals. We treat each spherical harmonic of the noise as 
an independent Gaussian random variable with variance 
described by the observed power spectrum [T3]. When 
data, such as the observed velocities V obs , are the sum 
of a signal V BH and a noise V n , the likelihood that the 
signal is present in V obs is the same as the likelihood that 
V — V BH is a Gaussian realization of the noise in the 
absence of a signal [19]. This allows us to determine the 
signal-to-noise ratio S/N for a given event: 
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Contributions to the S/N are greatest at harmonic de- 
grees / and frequencies v where the ratio of the numerator 
in the summation (left panel of Fig. [2]) to the denomina- 
tor (right panel of Fig. [2| is maximized. 

A PBH on a radial orbit freely falling from infinity 
(abbreviated as RFF hereafter) excites solar oscillations 
with S/N = 240 for m BH = 10~ 10 M Q . This is a conser- 
vative estimate for several reasons. Our simulations last 
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FIG. 4: Signal-to-noise ratio S/N as a function of maximum 
harmonic degree / for 1O- 1O M PBHs on the orbits shown in 
Fig. [3] The total S/N for each orbit is listed on each curve. 



only for t sim — 8 h and thus neglect contributions to the 
signal from t > t S i m . However, most of the signal comes 
from modes with frequencies v above the acoustic cutoff 
frequency p3] . These modes are only partially trapped in 
the solar interior, escaping into the corona on the sound- 
crossing time t ~ 2 J dR/c s = 1.96 h, where c s (R) is 
the sound speed as a function of solar radius. Simula- 
tions with t s i m > r should therefore capture most of the 
signal. Excluding the core (R < 0.24 Rq) from our simu- 
lations further reduces the signal, by both neglecting the 
energy deposited in this region and absorbing waves that 
reach the inner boundary. This absorbing inner bound- 
ary condition removes modes with v > 0.2Z mHz, as can 
be seen by the dearth of power in this region in the left 
panel of Fig. [2] The loss of these modes is significant be- 
cause the noise has very little power in this region as seen 
in the right panel of Fig. [2] Preliminary simulations with 
a core radius of 0.1 Rq suggest that the true S/N could 
be greater by a factor of 2 or more. These arguments 
imply that S/N = 240 is indeed a conservative estimate 
for m BH = 1O- 1O M = 2 x 10 23 g. Since S/N cx m BH , 
the minimum detectable PBH (S/N ~ 1) on the RFF 
orbit will have a mass m BB — 10 21 g. 

We have performed simulations with PBHs on three 
additional orbits: radial infall with a typical galactocen- 
tric velocity of 220 km/s at infinity (RGC), a parabolic 
orbit that maximizes the time spent at R < Rq (MT), 
and a parabolic orbit that barely skims the solar surface 
(SS). These orbits are shown in Fig. [3] The S/N for 
all four simulations as a function of the maximum / in- 
cluded in the summation in Eq. ([2| is shown in Fig. [4] 
Most of the power comes from large scales (/ < 100), since 
modes with higher I are evanescent rather than oscilla- 
tory deep in the solar interior where much of the energy 
is deposited [T3]. The small difference between the RFF 
and RGC simulations, both of which have no angular mo- 
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mentum (L = 0), suggests that S/N depends weakly on 
the orbital energy E. The RFF and MT orbits, both of 
which are parabolic (E = 0), demonstrate that the S/N 
depends more strongly on L. Our preliminary simula- 
tions with a core radius 0.1 Rq suggest an even stronger 
dependence on L, since the S/N of the RFF simulation 
doubles while that of the MT simulation remains nearly 
unchanged. The surface-skimming (SS) orbit has very 
low S/N, implying that the PBH must penetrate deeply 
into the solar interior to excite appreciable oscillations. 

One of our primary motivations to search for PBHs is 
the possibility that they constitute the cold dark mat- 
ter required by cosmology. The Milky Way resides in 
a dark-matter halo whose local density is approximately 
Pdm — 5 x 10 -25 g/cm 3 in the solar neighborhood. The 
Sun orbits the Galactic center with a velocity Vq ~ 220 
km/s. If PBHs with mass uibh and velocity dispersion 
a = vq/\/2 constitute the dark matter, the differential 
rate at which PBHs with specific energy E and angular 
momentum L encounter the Sun is [20] 

jm_ = 2p DM L [¥ c _ (2E+v%)/2(T z ginh / v Q V2E \ _ 
dEdL m B HVQ(T \ E \ a 2 J 

(3) 

PBHs with L < L max = ^2Rq{GM + ER ) have peri- 
centric distances less than Rq, implying that the total 
rate at which PBHs pass through the Sun is 
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If all PBHs passing through the Sun were detectable, 



Eq. Q would provide the event rate of our proposed sig- 
nal. Since the S/N depends on msn, E, and L, the true 
event rate can be found by weighting the integrand with 
a Heaviside step function that vanishes when the S/N is 
below a chosen threshold. Unless there is a considerable 
enhancement in the local PBH density beyond that ex- 
pected for Galactic dark matter, PBHs are unlikely to be 
discovered by solar observations alone. 

Fortunately, asteroseismologists study oscillations in 
stars other than our Sun. Such oscillations have been 
observed by the CoRoT (Convection Rotation and Plan- 
etary Transits) [5T] and Kepler satellites [22] . Resolution 
limits these instruments to disk-averaged observations 
sensitive to only the lowest I. Although this reduces the 
S/N of each event, the large number of stars that can be 
continuously monitored could greatly increase the total 
rate of detectable events. Future missions like the pro- 
posed Stellar Imager [23J may even resolve stellar disks, 
allowing S/N comparable to that for solar events. Fur- 
ther work is needed to establish that PBHs can excite 
detectable oscillations in stars with structures very dif- 
ferent from our Sun. If particle dark matter is not de- 
tected directly or discovered at the Large Hadron Col- 
lider, searches for alternative candidates like PBHs must 
be considered. We believe that asteroseismology might 
play an important role in these efforts. 
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